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Abstract 

We investigate the equilibrium morphology of a finite volume of liquid placed on two parallel rigid fibers of 
different radii. As observed for identical radii fibers, the liquid is either in a column morphology or adopts a drop 
shape depending on the inter-fiber distance. However the cross-sectional area and the critical inter-fiber distance 
at which the transition occurs are both modified by the polydispersity of the fibers. Using energy considerations, 
we analytically predict the critical inter-fiber distance corresponding to the transition between the column and the 
drop morphologies occurs. This distance depends both on the radii of the fibers and on the contact angle of the 
liquid. We perform experiments using a perfectly wetting liquid on two parallel nylon fibers: the results are in 
good agreement with our analytical model. The morphology of the capillary bridges between fibers of different 
radii is relevant to the modeling of large arrays of polydisperse fibers. 


1 Introduction 

The wetting of rigid and elastic fibers is present in 
man y natu ral and engineered systems: paper S and tex¬ 
tiles, EEHl fi lters , fog -harvesting nets, ^ human and 
animal hair or, more recently, microfabri- 

cated systems. Indeed, various biological or tech¬ 

nological materials are modeled as a network of fibers 
that are parallel or rando mly oriented. The presen ce 
of ambient humidity, the impact of a drop EHU or 
the addition of a fluid leads to the formation of capillary 
bridges similar to to t hose observed for liquids in gran¬ 
ular materials. The ability to model networks of 

fibers partially saturated with a liquid is crucial to deter¬ 
mine their mechanical properties. A key point for such 
modeling is the ability to describe the liquid morphology 
between fibers. 

Fundamental studies have considered the wetting ESEU 
and t he equi librium shape of a drop lying on a single 
fiber. Based on energetic considerations, it has 

been shown that a drop of liquid deposited on a fiber can 
either be in an axisymmetric barrel morphology engulfing 
the fiber or an asymmetric clamshell morphology with 
the droplet sitting on the side of the fiber. ^ Methods 
involving electrowetting effects have also been considered 
to tune the shape of the drop. EHEU 

However, understanding the wetting of complex fi¬ 


brous media first requires modeiing the morphoiogy of 
capiiiary bridges between severai fibers. Princen has 
characterized the coiumn morphoiogy that the iiquid 
adopts on two or more paraiiei fibers of the same radii. ^ 
In the coiumn morphoiogy, the iiquid is confined between 
the two fibers and has a constant cross-section aii aiong 
the coiumn. If some liquid is added or removed from 
the column, it will grow or shrink in length while its 
cross-section is not modified. Princen provided an an¬ 
alytical description of the shape of the liquid cross sec¬ 
tion in the column morphology, which makes amenable 
a prediction of the existence of a liquid column at small 
inter-fiber distances. Above a critical distance between 
the fibers, the liquid collects in a drop, which ove rspills 
the two fibers, as can be observed experimentally. EHHl 
The experiments with a drop on two rigid parallel fibers 
have shown that the transition between the drop and 
column morphologies is hysteretic. More recently, the 
equilibrium morphologies of the liquid on crossed fibers 
have also been studied both analytically and experimen¬ 
tally. ^ In this situation, a third morphology is possi¬ 
ble: the mixed morphology where a drop on one side 
together with a column of liquid on the other side co¬ 
exist. Whereas these different studies considered one or 
two rigid fibers, the competition between the elasticity 
of the fibers and the capillary effects has also been in¬ 
vestigated, ^ showing that a bundle of fibers organizes 
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into clusters. I ^ 5 | 36 | 37 | 

To the best of our knowledge, all of these studies de¬ 
voted to modeling the wetting of fibers focused on a pair 
or a large assembly of identical fibers. Yet, many biolog¬ 
ical or industrial systems, such as glass wool, involve a 
network of polydisperse fibers, i.e., fibers with different 
radii. In previous work, the capillary bridge between two 
or more spheres of different radii has been described, as 
well as its influence on the volume or the force exerted 
between the spheres but the particular situation 

of fibrous material remains to be considered. Indeed, the 
difference in fiber size is expected to affect the morphol¬ 
ogy and stability of the capillary bridge. 

In this paper, we build on the seminal study conducted 
by Princen ^ and we consider the morphology of a small 
volume of liquid placed on two parallel fibers of different 
radii, ai > a2. In principle, this study could be per¬ 
formed for fibers randomly oriented but the calculation 
would be more cumbersome without adding to the phys¬ 
ical picture. As a result, we focus on two parallel fibers 
with different radii and present an analytical description 
of the shape of the liquid cross-section in the column 
morphology. We also determine the critical inter-fiber 
distance at which the column to drop transition occurs. 
The results of the analytical model are then compared 
with experimental data obtained with a perfectly wet¬ 
ting fluid. The good agreement between the predictions 
and experiments validates our approach and our analyt¬ 
ical modeling. 


2 Analytical modeling 

We consider two rigid fibers of radii ai and a2 (ai > a2) 
and constant cross-section whose respective axes are de¬ 
noted {Oi z) and {O2Z). The fibers are parallel and 
separated by a distance 2d (see figure [^a)). If the dis¬ 
tance between the fibers is smaller than a critical distance 
2dc, we assume that a wetting liquid, with contact angle 
Oe < 90 °, deposited on the fibers spreads into a column 
shape. As shown by Princen,!^ the column has a con¬ 
stant cross-section except at the two terminal menisci 
whose analytical description remains difficult. However, 
for a sufficiently long column, assuming a constant cross- 
section provides a good approximation. The perpendic¬ 
ular cross-section of a liquid bridge between two fibers 
of radii ai and a2 {ai > a2) is schematically represented 
in figure [^b). The shape of the cross-section is charac¬ 
terized by its radius of curvature IZ and the wrapping 
angles of the fluid around the large and small fibers, ai 
and a2 respectively. 
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Figure 1 : (a) Top view schematic of the system composed 
of two fibers of radii ai and a2 separated by a distance 
2d. (b) Cross-section of a liquid column where ai and 
a2 denote the wrapping angles of the fluid around the 
fibers and is the liquid-fiber contact angle. 


ai and 0^2 is described by 

sin ai _ O2 C _ a2 + CA2 . . 

sin 0^2 OiC ai^CAi 

In addition, the distance between the axes of the two 
fibers Oi O2 is given by: 

O 1 O 2 = ai+a2 + 2d = P0i-|-P02 

= (ai + CAi) cosQ^i + (a2 + CA 2 ) cos(a2, 

( 2 ) 


which leads to a relation between o^i, 0^2 and d: 

CAi cosai+CA2 cos 0^2 = ai(l—cos(ai)+a2(l—cos a2)-\-2d. 

( 3 ) 

The triangle A1O3A2 is isosceles with IZ = O3A1 = 
O3A2. In addition, O3C defines the bisector of angle 
A1O3A2. As a result, CAi and CA2 are equal: 

ai(l - cosQ^i) + a2(l - cos(a2) + 2(i 

GAi — CA2 — -. 

cos ai + cos 0^2 

( 4 ) 

The radius of curvature 7Z of the cross-section is given 
by 

7 ^ = O3A1 = CAi ( 5 ) 

smp 

where 77 and P are two angles defined in figure Ge¬ 
ometrically, we find that 77 = {t: — ai — a2)/2 and 


2.1 Geometrical determination 

Because two inscribed triangles share a common side CP 
(figure Bb)), the relation between the wrapping angles 
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Figure 2: Focus on the cross-section of a liquid column. 

(3 = (tt — 0^1 — 0^2 — 2 ^£;)/ 2 . Using the radius of the 
large fiber, ai, as the characteristic length scale and the 
ratio of the fiber radii, denoted F = a2/ai, equation ^ 
can be written: 

^ IZ (1 — cos (^i) + F (1 — cos 0^2) + 2 d 
ai cos((ai + Oe) + cos((a2 + Oe) 


where we assume that the increase in volume leads to 
an increase of length of the column dV = ^dL. In this 
relation, 7 is the liquid/vapor surface tension, P and Pq 
are the pressure values inside and outside the column, 
respectively. In addition, the increase in length of a col¬ 
umn results in the creation of liquid-fiber and liquid-air 
area while reducing the air-fiber area. This process is 
associated with an energy variation 

dE = [{AiBi + A2B2) {'jsL — 7 sv) + (^1^2 + B1B2) 7] dL, 

(11) 

where ^sl and jsv are the solid/liquid and solid/vapor 
surface tensions, respectively. In this relation AiBi + 
A2B2 = 2 aiai + 2a2a2 denotes the region of the fibers 
initially non-wetted by the liquid and ^1^42 P B1B2 = 

{ 7 T — 20 E — oti—a2)Tl describes the creation of a liquid- 
vapor interface. Using the Young’s relation, we have 
ISL — Isv = 7 cosOe- These results lead to an energy 
change 


where d = d/ai. We also need to determine the cross- 
sectional area A of the liquid column. A calculation 
based on the geometry of the system is detailed in the 
appendix and leads to 

A 

—2 = 'JZ[cos{ai P Oe) P cos{a2 P Oe)] (sinai + F sin 0 ^ 2 ) 

CLi 

— {ai — sin 0^1 coso^i) — F^ (0^2 — sin 0^2 cos 0^2) 
OL\ -|- 0^2 + 2 Oe \ f OL\ + 0^2 T 2 Oe 


dE = 2 (7r — 2 ^£; —0^1—0^2) 7^—2 ((aiaiT(a2a2) cos ydL. 

( 12 ) 

The first term in this expression is proportional to the 
capillary force at the air-liquid interface and the second 
term is proportional to the force exerted by the fiber on 
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Finally, a last relation between ai and 0^2 can be ob¬ 
tained: 


{d P 1 ) tan = {dpr) tan 


(8) 


d = 


F tan((a2/2) — tan((ai/2) 
tan(ai/2) — tan((a2/2) 


( 9 ) 


the liquid. Balancing the expressions ( 10 ) and ( 12 ) for 
the added volume dV and using the expression of the 
cross-sectional area, we obtain the equation: 


211 - 


2 r TT 0^1 H- 0^2 + 2 Oe 


, 'ai P a2p20E 

P sm I -^- I cos 


L2 2 

di P (12 P ‘2‘ Oe 


pn 


Alternatively, this expression can also been rearranged 
to obtain d(r, q^i, 0^2): 


2 J \ 2 

(sinofi + F sin0^2) (cos((ai + Oe) P cos((a2 + Oe)) 

—2 {ai + F0^2) — {di — sin0^1 coso^i) 
—F^ {a2 — sin0^2 cos 0^2) = 0 


( 13 ) 


In the following, we focus on a perfectly wetting fluid, 
i.e. with a contact angle 6>e = 0 °. In this situation, the 
previous expression becomes 


2.2 Energy balance 

Following Princen’s approach,!^ we neglect the contri¬ 
bution of the radius of curvature in the direction parallel 
to the fibers. Indeed, far from the menisci, located at 
both ends of the column, the interface is nearly flat and 
its contribution is small with respect to the radius of 
curvature in the cross-section. The variation of energy 
associated with the addition of an infinitesimal volume 
dV to the column is equal to the work dIU done on the 
system as a result of a pressure difference with the at¬ 
mosphere: 

dE = dW = {P-PQ)dV = -lAdL, ( 10 ) 




TT Q^i + 0^2 . / 0^1 + 0^2 \ / Q^i + 0^2 

-h sm - cos - 

2 2 V 2 y V 2 


PP (sinofi+F sin0^2) (cos ai + cos 0^2) 
—2 (ofi + F 0^2) — {di — sin 0^1 cosq^i) 
—F^ {d2 — sin0^2 cos 0^2) = 0 . 


( 14 ) 


2.3 Numerical solution 
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Considering a perfectly wetting fluid {Oe = 0 °), the 
shape of the cross-section is fully characterized by o^i. 
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Figure 3 : Analytical shapes of the cross-section for dif¬ 
ferent radii ratios F = a^jdx and values of the wrapping 
angle around the largest fiber, ai. 


0^2, 'IZ and d for a given value of the fiber aspect ra¬ 
tio F = a^jdx. Typically, we want to vary ai G [ 0 , tt] 
and successively obtain the corresponding values of 0^2, 
d and IZ. The calculations are performed using a custom- 
written Matlab cod^ The wrapping angles ai and 0^2 
are increasing when increasing the inter-fiber distance d. 
When the liquid wraps around the fiber, ai increases 
until it reaches a maximum ai = tt. When = tt, the 
liquid can not wrap around the fiber anymore therefore 
leading to a transition from a column morphology to a 
drop morphol ogy, s imilarly to the situation of identical 
parallel fibers. ESEU addition, equation ^ shows that 
for ai = TT, we have 0^2 = tt. Substituting these values in 
equation ( 14 ) and ^ we obtain: 


4 = G = 

ai 


( 15 ) 


Considering two identical fibers leads to F = 1 and 
dc = a/2 , which recovers the result previously obtained 
for identical fibers. ^ 


Using the solution of equation ( 14 ), we can determine 
the shape of the cross-section for various radii ratios F 
and for increasing values of the wrapping angle ai. In 
figure we have represented three situations: F = 1, 
i.e. identical fibers, F = 1/2 and F = 1 / 10 . We observe 
that 0^2 increases faster than ai for F < 1 and that the 
column can become convex outward for ai < 90 °. In 


^Available in Supplementary Material (ESI) 



Figure 4 : Evolution of L/V = 1/.4 with the inter-fiber 
distance 2 d for various fiber aspect ratios F = 1 (black) 
and F = 0.5 (red), F = 0.1 (green) and F = 0.05 (blue). 
The circles at the end of each curve and the sketches 
of the cross-sectional area correspond to the maximum 
inter-fiber distance for the aspect ratio F considered. 


addition, the shape of the cross-section associated with 
the maximum area A is obtained for o^i = 0^2 for any 
value of F; only the evolution of the shape with ai is 
modified with the radii ratio. 

For a given value of F and an inter-fiber distance d, we 
can determine the wrapping angles {ai and 0^2) and 1 Z, 
which leads to a unique value of the cross-sectional area 
A. Consequently, we define the length of a column L = 
L/di for a given volume of liquid V = V/di^ deposited 
on the fiber. Following Protiere et al.,l^ we report the 
theoretical rescaled parameter L/V = 1/^ in figure]^ 
The influence of the radius ratio F can be observed both 
in the value of the critical inter-fiber distance dc and in 
the value of 1 /^. The evolution of L/V is non trivial: for 
a fixed d, the values of and 0^2 lead to a smaller cross- 
sectional area for large radius ratio therefore leading to 
longer column. 

In the next section, we compare our analytical results 
to experiments performed with different radii ratios F 
and a perfectly wetting fluid. 


3 Experiments 

A typical experiment is performed by depositing a small 
amount of silicone oil (viscosity 5 cSt), which is a per¬ 
fectly wetting fluid (d^; = 0), on a pair of nylon fibers 
of circular cross-section and respective radii (ai;a2) G 
[ 50 , 225 ] /im. Each fiber is clamped on two horizontal mi¬ 
cro controllers (PTl, Thorlabs) that allows very accurate 
increases of d by step of 5 jam. To minimize the effects of 
gravity, we deposit a small amount of liquid V G [ 1 , 6] /i^ 
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(b) 

Figure 5 : Evolution of a drop of 5 cSt silicon oil placed 
on two parallel nylon fibers as the distance d between the 
fibers is increased for (a) equal radii ai = a2 = 225 /iui 
(V = 6 lit) and (b) different radii ai = 225 /im and a2 = 
75 yum (V = 3 lit). The dashed line indicates the column 
to drop transition. Scale bar is 5 cm. 


Figure 6: Evolution of the rescaled length L = Ljai of 
5 cSt silicon oil placed on two parallel nylon fibers as a 
function of the rescaled distance d = d/ai between the 
fibers for F = 2/3 (ai = 225 /im) and for three differ¬ 
ent rescaled volumes V = 88, V = 263 and V = 527 . 
The solid line is the prediction of the model described in 
section and the circles are the experimental measure¬ 
ments. The insets show the transition between drop and 
column for V = 263 and V = 527 (top pictures: top view; 
bottom pictures: side view). 


using a micropipette and we neglect the hysteresis in¬ 
duced by the gravitational effects that deform the drop 
into a hemispherical shape. ^ Indeed, a deviation from 
the analytical prediction is observed when the drop to 
column transition is reached by reducing the inter-fiber 
distance 2 d. ^ We therefore start from a small inter-fiber 
distance 2 d and increase this distance using the microm¬ 
eter stage. Pictures are taken from both the top and side 
views to measure the length of the column and to discrim¬ 
inate the column and drop morphologies. In particular, 
in the drop morphology, the width of the drop is larger 
than the spacing of the fibers, which implies that a bead 
of liquid is located outside the fibers. Through direct 
measurement, we ensure that the distance between the 
fibers increases in steps of 10 /im. After increasing the 
inter-fiber distance, we typically wait for a few minutes 
so that the equilibrium state is reached prior to imaging. 

Typical visualizations for a pair of identical fibers and 
a pair of different fiber radii are reported in figure |^a,b), 
respectively. We observe that the wetting length L de¬ 
creases in both situations when increasing the inter-fiber 
distance 2 d. However, the critical distance at which the 
column to drop transition occurs, dc, is different for the 
two cases, as predicted by the analytical model. We also 
observe that the menisci on both extremities of the col¬ 
umn become non-symmetric with fibers of different radii 
(figure l^b)). 

The values of the wetting lengths L measured in both 


situations are compared to the analytical predictions de¬ 
rived in the previous section. In figure we report the 
evolution of L for various rescaled volumes and non¬ 
similar fibers. The analytical model captures the evo¬ 
lution of the length, which further confirms our assump¬ 
tion to neglect the menisci at the ends of the column. 
For large volumes, the column to drop transition leads 
to a jump in the measured length, whereas for small liq¬ 
uid volumes this transition remains smooth. Therefore, 
to estimate the critical distance dc at which the column 
to drop transition occurs, we rely on measurements per¬ 
formed with sufficiently large volumes. We should em¬ 
phasize that, for large volumes, a hysteresis between the 
column to drop and drop to column transition is observed 
similar to the report of Protiere et al. ^ To predict an¬ 
alytically the value of the hysteresis would require con¬ 
sideration of the energy in the column morphology and 
in the drop morphology. However the analytical predic¬ 
tion of the drop morphology remains complicated and 
would necessitate performing numerical simulations that 
minimize the energy in th e presence of gravity using, for 
instance. Surface Evolver. Pi | 42 | 43 | 44 ] 

Starting from a column shape, we can measure the 
critical inter-fiber distance 2 dc at which the column to 
drop transition occurs and compare these experimental 
results to our analytical prediction. Analytically, we find 
that for a perfectly wetting fluid the transition occurs 
for dc = a/ 2T (equation [l^. We have performed experi- 
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Figure 7 : Evolution of the critical distance for the column 
to drop transition when the inter-fiber distance 2 d is in¬ 
creased for various pair of fibers of radii ai and a2. The 
squares are experiments performed with ai = 225 /im 
and: a2 = 50 /im and V = 1 jj^i (blue), a2 = 75 /im and 

V = 1 jj^i (green), a2 = 100/im and V = 3 fii (brown), 
a2 = 150 jam and V = 3 jai (magenta), a2 = 225 jam and 

V = 6 jii (cyan). Hollow circles shows experiments with 
with ai = 150 jam and: a2 = 50 jam and V = 1 jai (blue), 
a2 = 75 jam and V = 1 jai (red). The purple star is data 
obtained from Protiere et al. for ai = a2 = 125 jam and 

V = 4 jai (blue). ^ The continuous line is the analytical 
prediction ([T^. 


ments with pairs of fibers ai and a2 including the situa¬ 
tion described by Princen, ai = a2 leading to P = 1 . 
A summary of our experimental results is presented in 
figure We observe that for all values of F the analyti¬ 
cal prediction captures well the morphological transition, 
which further confirms our analytical modeling. 


4 Discussion and conclusion 


the distance y to the “crossing” point: E 3 


d{y) 


d{y) 

ai 



(i + r)2 (1 + r) 

4 2 ’ 

( 16 ) 


where S is the tilt angle between the fibers. It follows 
that the area of the cross-section depends on the varying 
inter-fiber distance 2 d(y) according to the relation ([t]). 
Therefore, we can determine A for all values of d < \/^. 
The condition dc = V^ 2 T predicts the maximum length 
of the column Lc- The length of the column for a given 
volume of fluid V deposited on the crossed-fibers satisfies: 


V = J ^A{y)dy. ( 17 ) 


If L > Lc the fluid is either in a mixed morphology or 
in a drop shape whereas for L < Lc the fluid is in a 
column morphology. With this geometric picture that 
characterizes the morphology, we could also determine 
the liquid-air area and thus the ev olutio n of the drying 
rate for the different morphologies. ESEH 

Similarly, the calculation could be extended easily to 
fibers of identical radius but with different surface prop¬ 
erties, i.e. if the contact angle of the fluid is different for 
the two fibers 9 e,i ^ 0 e, 2 - Similar behavior to the results 
presented in this paper is expected to happen: namely, 
the wrapping angle on one fiber, o^i, is different from the 
wrapping angle on the second fiber 0^2. It will modify the 
transition between the drop and column morphologies as 
well as the critical inter-distance dc compared to fibers 
of identical radius and the same surface properties. 

To conclude, we have taken a step toward the de¬ 
scription of wetting in array of fibers by considering the 
liquid morphology for fibers of different sizes. Work is 
in progress to better describe the distribution of liq¬ 
uid in a large network of fibers. First, as emphasized 
above, numerical simulations, for instance using Surface 
Evolverl^, are needed to predict the shape of a drop 
with and without gravity and determine the associated 
energy. Einally, more complex interactions, i.e., between 
three or more fibers, should be considered. 


In this paper, we have presented an analytical model, 
confirmed by our experiments, of the column morphology 
and of the critical inter-fiber distance for the column to 
drop transition. These results extend a previous analyt¬ 
ical study devoted to identical fibers to the situation 
of fibers of different radii. 

Although the present paper focuses on the situation of 
parallel fibers, since it constitutes a simpl e case that has 
been well studied in the literature,®^ the analytical 
calculation is straightforward to extend to the situation 
of crossed-fibers. Indeed, for crossed fibers, the main 
difference is that the inter-fiber distance 2 d depends on 
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A Appendix A: cross-sectional 
area of the liquid column 


To determine the area of a cross section ABCD, we con¬ 
sider the area of the trapezoid shown in Fig. 


^ _{AB^ CD) A 

— 7 ^ 


(18) 


Substituting this expression in relation ( 25 ) leads to 


A 

—^ = 7 Z[cos{ai + Oe) + cos(a2 + Oe)] (sinai + F sin0^2) 
— {ai — sin 0^1 cosai) — F^ (0^2 — sin 0^2 cos 0^2) 


-2 7 ^^ 


TT 

-Sin 

L2 


cti + ct2 + 2 Oe 


cos 


cti + 0^2 + 2 Oe 


cti + ct2 + 2 Oe 


with AB = 2 ai sin^i, CD = 2 a 2 sin 0^2 and A = 2 d-\- 
ai (1 —cosai)+a2 (1 —cos 0^2). We then remove the areas 
SI, S2, S 3 and S 4 , which are given respectively by: 


Considering two identical fibers of radii ai = a2 = a 
leads to ai = 0^2 and the cross-sectional area simplifies 
to 


^2 

^3 


^4 


cil [<^i — sino^i cosQ^i] 

^ [0^2 — sin 0^2 cos 0^2] 

’TT CXi 0.2 2 0E 

.2 2 


Ci 2 


— sm 


^3. 


Oi 02 2 0E 


cos 


( 19 ) 

( 20 ) 


+ 0^2 + 2 

( 22 ) 


2 ai 



A 

= AIZ Ano cos{o + 6>£;) — 2 {o — sino cos o) 


-2 7 ^^ 


--o - Oe - An {o A Oe) cos {o A Oe) , ( 26 ) 


which is the expression established by Princen for a pair 
of identical fibers. ^ 
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Eigure 8: Schematic of the cross-section to determine its 
area. 


Therefore, we obtain the cross-sectional area 


A = (ai sin0^1 + a2 sin0^2) [ai (1 — coso^i) + a2 (1 — cos0^2) + 2 d] 
—ai^ {oi — sin0^1 coso^i) — a2^ {02 — sin0^2 cos0^2) 


-2 7 ^" 


TT 

L 2 ~ 


oi A 02 A 2 Oe 


cos 


0^1 + 0^2 + 2 


0^1 + 0^2 T 2 Oe' 


( 23 ) 


In addition, relation ^ can be rewritten 

(1 — cosQ^i) + r(l — cos 0^2) + 2 d = lZ[cos{oi + Oe) + cos((a2 + Oe)]- 

( 24 ) 


(25) 
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